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ABSTRACT 

We present an analytical study of the statistical properties of integrated emission and velocity centroids for a 
slightly compressible turbulent slab model, to retrieve the underlying statistics of three-dimensional density and 
velocity fluctuations. Under the assumptions that the density and velocity fields are homogeneous and isotropic, 
we derive the expressions of the antenna temperature for an optically thin spectral line observation, and of its 
successive moments with respect to the line of sight velocity component, focusing on the zeroth (intensity or 
integrated emission /) and first (non-normalized velocity centroid C) moments. The ratio of the latter to the 
former is the normalized centroid Co, whose expression can be linearized for small density fluctuations. To 
describe the statistics of /, C and Co, we derive expansions of their autocorrelation functions in powers of den- 
sity fluctuations and perform a lowest-order real-space calculation of their scaling behaviour, assuming that the 
density and velocity fields are fractional Br ownian motions. We hence confi rm, within the scope of this study, 
the property recently found numerically bv iMiville-Deschenes et ail l|2003a|) that the spectral index of the nor- 
malized centroid is equal to that of the full velocity field. However, it is also argued that, in order to retrieve the 
velocity statistics, normalization of centroids may actually not be the best way to remove the influence of den- 
sity fl uctuations. In this respect, we discuss the modified velocity centroids introduced hv lW-arian fc Esnuivell 
l|2003h as a possible alternative. In a following paper, we shall present numerical studies aimed at assessing the 
validity domain of these results. Appendices lA"l to IT)1 are only available online at EdP Sciences. 

Subject headings: ISM: structure - Methods: analytical - Turbulence 

1 Introduction 



The proper exploitation of astronomical observations requires one to deal with several problems to accurately 
describe the objects and processes under study. In particular, regarding the physics of the interstellar medium 
(ISM), it should be stressed that spectral emission data depends on the velocity field solely via its component 
along the line of sight, through Doppler shifts. Furthermore, this information is necessarily integrated along the 
line of sight, and radiative transfer leads to expression s in w hich the contributions of the density and velocity 
fields are mixed in a complex way l|Hegmann fc Kegel feoom. Hence, to describe the physical conditions and 
processes in the ISM, one has to rely on a single number (e.g. antenna temperature) for any given direction 
in the plane of the sky and any given velocity along the line of sight. Although the comparison of antenna 
temperatures for various tracers helps, it is necessary to solve an inverse problem to have access to the three- 
dimensional properties of the medium, such as density and velocity, and compare them with various models, 
for instance to assess the roles of gravity, magnetic field or turbulence. 



1 



Indeed, it is now recognized that the different components of the ISM are subject ed to turbulent motion. Thi s has 
been observed in the ionized gas llVan Langevelde et allll992h. in Hn regions llO'Dell fc Ca,sta,nedaLlT987h and 
in the neutral atomic phase l|Spicker lTF^rtzh ige^{jl98q; [Miville-peschenes et a^|2003 bl). but molecular cloud 
studi es are by far the most numerous (see e.g. lKleiner fc Dickmanlll985allKitamura et al I ll993UMiesch fc Ballvl 
Il994h . Estimates of Reynolds numbers from molecular viscosity in these clouds are of the order of 10 8 , consistent 
with turbulent flows l)Chandrasekharlll949(l . Moreover, molecular lines in this phase exhibit suprathermal widths 
over a wide range of spatial scales, from a few km s _1 in small dark clouds to a few tens of km s _1 in giant 
complexes, while the thermal dispersion is only of a bout 0.3 kins" 1 for molecular hydrogen at T = 10 K. These 
linewidths scale as a power-law of the cloud's si ze llLarsorill98l . with an exponent close to that predicted by 
the classical theory of turbulence l|Kolmogorovl Il94ll) . Such r andom motio ns within molecular clouds may in 



turn account for a number of other properties of spectral lines feakeit lT97fn. as well as provide support against 
gravi tational colla pse, explaining the fact that the lifetime of molecular clouds is larger than their free-fall 
time (|Scalol Il985|) and t hat the star formation rate is therefore much smaller than predicted by gravitationally 
collapsing cloud models l)Zuckerman fc Evanslll974j) . 

Because of this interplay between random motion and many physical processes at work in the ISM in general 
an d in the molecular phase in part icular, one needs to accurately describe these turbulent flows (see the review 
bv lVazauez-Semadeni et aill2000() . As noted earlier, this has to be done using the antenna temperature which 
represents the emission from a given direction and at a given velocity, and so a number of methods were devised 
to derive the statistical properties of the thre e-dimensional fields from tho se of the observational data. Among 
these, the velocity channel analysis (VCA) ofL azarian fc PogosvarJ iTjlKHH is based on an analytical derivation 
of the properties of channel maps with varying velocity widths. It may however prove difficult to apply to 
ac tual observations, as show n bv IMiville-Deschenes et aTl l|2003a[) . The modified velocity centroids (MVC) 
of lLazarian fc Esauivell l)2003[) are a recent promising attempt to reduce the influence of density fluctuations in 
the statistical properties of centroids, although it can be argued that they are only defined through their structure 
function. As a final example of velocity statistics retrieval methods, principal component analysis (PCA), which 
works on the full position-position-velocity cubes, is meant to decompose data onto an orthogona l basis and 
derive properties of the velocity field at each scale, as calibrated numerically bv lBrunt et alJ l)2003[) . The main 
objective of these works is to relate the scaling behav iour observed i n the two-dimensional maps to scaling laws 
inferred for the three-dimensional fields. For instance. IStutzki et alJ l)l998|) showed that for optically thin media, 
the spectral index of the integrated emission map is the same as that of the ful l density field, provided t hat the 
depth probed is larger than the transverse scales considered. In a recent work. IMiville-Deschenes et"al . (2003a) 
used numerical simulations to show that the same is true for the normalized velocity centroid with respect to 
the three-dimensional velocity field. However, this latter result lacks theoretical support, and it is therefore 
the goal of this paper to present an analytical study aimed at clarifying the relationship between the velocity 
centroids and the velocity field, within a simple turbulent cloud model. 

Given the observational data, one may derive moments of the antenna temperature profiles, each of these 
moments yielding a potentially informative two-dimensional map. For instance, the zeroth moment is the inte- 
grated e mission, or intensity, while the first moment is the velocity centroid, which ca n also be normalized to the 
intensity ijMiinchl . Il958t iDickman fc Kleinerlll985bl:lMiville-Deschenes etall l2003a|) . For an optically thin line 
and uniform excitation conditions, the non-normalized centroid can be related to the cloud's total momentum, 
while the normalized centroid is a measure of av erage velocity within the m edium. While density fluctuations 
may bias the description of the turbulent motion l Lazarian fc Esauivell I2003T) . it is however commonly believed, 
and intuitively plau sible, that their effects are somewhat compensated by normalization. To properly assess 
these , and following TPickman fc Kleinerl l)l985b|) (see also lScalcl Il984t Ixitamura et all Il993t FMiesch fc Ballvl 
Il994h . we shall use autocorrelation functions of the moment maps and relate them to correlation functions of 
the underlying three-dimensional fields. To this end, we first describe the model we shall use and introduce the 
notations and assumptions in sectional We then present a brief summary of how moments of the line profile can 
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be related to the density and velocity fields within the medium (section |3J) . Section 0] contains a general study 
of the statistical properties of intensity, normalized and non-normalized velocity centroid maps as functions of 
the three-dimensional density and velocity fields' statistics. The equations obtained in the lowest order are then 
applied to the test case of fractional Brownian motion (fBm) density and velocity fields (section EJ. Section H3 
presents a discussion of the various results obtained with respect to earlier works. Our concluding remarks are 
given in section [7] and details on the calculations can be found in the appendix. 




1 1 1 1 

-D/2 D/2 z 

Figure 1: Notations used in the paper (left) and schematics of the turbulent slab model (right). A lowercase 
boldface letter stands for a three-dimensional vector, while the corresponding uppercase represents its projection 
on the plane of the sky (xOy). The slab is infinite in the x and y directions and is limited to Iq = [—D/2, D/2] 
in the z direction. The density p and line-of-sight velocity v distributions along a line of sight are shown on the 
right. The mean values po and vq are taken over the whole slab. 



2 The model 

Throughout the paper, boldface notations stand for vector quantities. A point in three-dimensional space is 
given by its position x, which can be written as x = (X, z) = X + ze z , where X is a two-dimensional vector 
in the plane of the sky and z marks the line-of-sight position of the point considered, e z being the unit vector 
along the line of sight. The three-dimensional separation between points Xi and X2 is written as r = X2 — xi and 
the separation between their respective lines of sight Xi and X2 is R = X2 — Xi. In short, three-dimensional 
vectors are written in lowercase, while vectors in the plane of the sky are written in uppercase. These notations 
are illustrated in Fig^ It should be emphasized that we shall only consider small scales on the plane of the 
sky, so that all lines of sight are parallel. 

Hereafter, our model is a slab of gas of width D, perpendicular to the line of sight, and of infinite transverse 
extensions 1 , with z = conventionally placed halfway through the slab (see Fig.QJ. Within the slab, the average 

1 This is in no way contradictory with the assumption that all lines of sight considered are parallel. The cloud is supposed to be 
infinite, but only a small fraction of it is observed and supposed to be representative of the whole. 
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density (that is, the average number of emitters per unit volume taken over the whole slab) is a constant noted 
Pq. As for the velocity v of the gas with respect to the observer, we shall write v to stand for its component 
along the line of sight, the mean value of which over the slab is vo- These averages (po and vq) are obviously 
constants, independent of the line of sight. To assume the most general case, we shall take vo to be nonzero 2 . 
The gas velocity within the rest frame of the cloud is simply <5v = v — voe z , assuming that there is no systematic 
transverse velocity. Furthermore, fluctuations of density along any given line of sight are supposed to be small 
compared to the average value po, in the sense that the variance of the density field should be less than p . It is 
also assumed that the turbulent flow w ithin the slab is h o moge neous and isotropic. These last hypotheses are 
to be understood in the strong sense of iMonin fc Yaeloml l|197q) : not only are the scalar field p and the vector 
field 5v homogeneous and isotropic, but the same should apply to the four-dimensional field (p, 5v). This will 
be useful in the derivation of our results in section 21 Outside the slab, the density and velocity are set to zero. 
We further assume that the observation is done in an optically thin transition, at a frequency for which the 
Rayleigh-Jeans approximation is valid, that uniform excitation conditions apply within the slab, and that no 
background radiation is present. 

With these assumptions in mind, the antenna temperature T a (X, u) representing the em i ssion fr om a given line 
of sight X at an observed velocity u can be written as the integral l(Dickman fc Kleineii ll985bTl . 

T a (X,it)= Jt cx (x)ko(x)(/)(v(x) - u)dz, (1) 

Id 

where Ijj is the segment [-D/2, D/2] over which the integration is performed. At each position x, T ex is the 
excitation temperature and kq is the integrated absorption coefficient. The normalized line profile function (f> 
is assumed to be symmetrical about zero and independent of x. For instance, if one only considers thermal 
broadening, <fi takes the form 



H w ) = -7= — ex P 

V ZTTC th 



2<4 



(2) 



with cth the thermal velocity dispersion. The shifted argument w(x) — u of in Eq. Q simply expresses the 
fact that the emission at position x is broadened around the local line of sight velocity f(x). Assuming uniform 
excitation conditions, T ex is a constant, and kq(x) is proportional to the local gas density p(x), so we can write 



T a (X,u) = J op(x)0(«(x) - u)dz, (3) 

In 

where a is a proportionality constant. The expression in Eq. shows that, even under simplifying assumptions, 
the task of extracting properties of the density and velocity fields from the observational data T a (X,it) is a 
difficult one. 



3 Velocity moments 

In order to obtain information about the velocity field properties from the data sets, one may compute the 
various moments of the antenna temperature profiles T a (X, u) for a given line of sight X. This is a logical step 
considering that knowledge of the moments of a random variable is equivalent to that of the full probability 
distribution. We shall therefore consider the moment maps W n (K) defined by 

W„(X) = Ju n T a (X,u)du. (4) 

2 The case vq = is only special in the way the calculations are performed, as the end result amounts to setting vo to zero in 
the general formulae. This is shown in appendix IbI 
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The integration is done over all velocities from — oo to oo, which poses no convergence problem since the line 
profile has a finite support 3 . Using Eq. J3J to express T a as a function of the local line profile and of the density, 



W„(X) = Ju n ajp(x)(f>(v(x) — u)dz du = ajp(x) J(v(x) — w) n 4>(w)dw 

I Id J Id I 



dz, 



(5) 



where we introduced the variable w = v(x) — u. Developing the integrand in the innermost integral, and since 
the local line profile cj> is assumed to be independent of the position x, we get the following expression 



W„(X) = c^(-l) fe C* 



fc=0 



p(x)v n - k (x)dz 



w 4>(w)dw 



= J2 c * hk 



fc=0 



P (x)v n - k (x)dz 



(6) 



where the C k are the binomial's coefficients and the hk are related to the moments of the normalized line profile 



h k = {-l) k ajw k (t>{w)dw. (7) 

It should be noted that the assumed evenness of 4> leads to ft-2p+i = 0. The first two moments arc of particuler 
interest, since, by definition, they are the total emergent intensity I and the non-normalized velocity centroid 
C, respectively, 



I(X) = W (X) = ajp(x)dz = aN(X) and C(X) = Wi(X) = ajp(x)v(x)dz. 

Id Id 



(8) 



Here, N(X) is the column density at position X. In the rest of this paper, we shall only consider the intensity 
and the velocity centroid. The use of higher order moments would theoretically give access to more information 
about the structures of the density and velocity fields, but, in practice, noise levels and limited spectral resolution 
may very well jeopardize their usefulness. 



4 Statistics of intensity and centroid maps 
4.1 Rationale of the computations 

We may use the line profile moments defined above to quantify the structure of the turbulent flow within 
the slab. Obviously, the zeroth moment Wo(X) = /(X) can only be used to derive statistics of the density 
field, as it is proportional to the column density N(X). The first moment or non-normalized velocity centroid 
Wi(X) = C(X) is the first appropriate quantity to study the velocity field, as can be seen from Eq. JSJ. 
However, C(X) represents the integration of a momentum, rather than of a velocity proper, and it appears 
that density fluctuations may affect the estimation of the velocity statistics from this map. To circumvent this 
problem, it is common to use normalized centroids Co(X), which are simply defined as the ratio C(X)/7(X). 
This is usually and empirically justified by the assumption that density fluctuations in C(X) are also present in 
I(X) and therefore somehow vanish from Co(X). In the simplified case where density does not fluctuate along 
the line of sight, such a reasoning is obviously correct, even if transverse density fluctuations are present, and 
normalization then only serves as a dimensionality factor. As far as we are aware, however, no analytical study 
exists of the influence of longitudinal fluctuations of density on velocity centroids, even in the quite simplified 
model used here. 

3 If any continuum emission is present, it should first be removed. 
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Now, in order to gain a better understanding of the underlying three-dimensional statistics of the velocity field 
through velocity centroids, it seems natural to consider the two-dimensional statistics of the centroid maps. 
Indeed, structural properties of the velocity and density fields should arise, under one form or another, in 
statistical measures performed on the maps C (X) and Co (X) . Of course, it is similarly reasonable to look for 
the influence of density structure in the statistics of the intensity map /(X). One such useful measure to be 
performed on the available two-dimensional maps is the autocorrelation function (ACF), which gives the mean 
degr ee of correlation between values of a field taken at points separated on the plane of the sky by a given vector 
R |Kleiner fc Dickmanl Il984|) . More precisely, the autocorrelation function Ayr of a field T is defined by 

4f(r) = <^( x )^( x + R)> . (9) 

where the brackets stand for a spatial average over position X in the plane of the sky. Hereafter, we shall 
consider the autocorrelation functions of the intensity and of both types of centroid, respectively noted Aj(R), 
Ac (R) and Ac a (R) , and compute them as functions of statistical measures of the three-dimensional density and 
velocity fields. In order to do so, we first separate mean and fluctuating contributions of density and velocity 
in the quantities involved, with p = po + Sp and v = vq + Sv. These expressions can be used to write the first 
two moment maps as 

/(X) = Jap(X,z)dz = ap D[l + y Px ] and C(X) = Jap(X, z)v(X, z)dz = ap v D [1 + y px + y Vx + y pV:K \ , 

(10) 



Id Id 

where we have introduced the following integrated fluctuation terms 4 



1 [ 6p(X,z) A 1 [MM, 1 f Sp(X,z) Sv(X^ 

D J po D J v D J po Vq 

Id Id Id 

The normalized centroid Co is then simply written as 

Co(x) -7(xy" Uo rr^ ■ (12) 

This last expression can be used to clarify the usefulness of the small fluctuations hypothesis. Indeed, for a 
perturbative method to be applicable in the computation of the autocorrelation function of the normalized 
velocity centroid map, one should be able to linearize Eq. I|12|l . and therefore we need to assume that \y Px \ < 1, 
so that the denominator can be expanded as a converging series, 



i-^x + y" x --- = E(-^x) n - (is) 



Such a condition is obviously achieved if local density fluctuations themselves are small in the sense, expressed 
in section |21 that the standard deviation of the density field a p should be smaller than the average density po, 
since 



\Vpy 



D J po 

Id 



\ 



1 f(5p{X,z)\\ _ I dp 2 



DJ\ po 

Id 




dz ^\Hr = ~- (14) 



4 In these integrated fluctuation terms, we chose to denote the position X on the plane of the sky as a subindex, rather than as 
an argument, to allow for concise equations in the subsequent developments. 



(i 



The identification of the mean squared density fluctuations along a single line of sight with the variance <7 p of 
the whole density field is based on the homogeneity hypothesis. Indeed, any line of sight should statistically 
represent the full field. Although the stronger assumption a p < po is not strictly necessary for the linearization 
mentioned earlier, it is useful in the expansion of the autocorrelation functions since, in this case, the hybrid 
integrated fluctuation term y pvx is of the same order in density as t/ px , 



\Vpvx\ < 



\ 



1_ j-^ SpjX, z) 



Id 



Pa 



dz 



1 f(5v(X,z) 



t'o 



Id 



dz ~ — — . 

Po v 



(15) 



where we used the Cauchy-Schwarz inequality and introduced the standard deviation a v of the line-of-sight 
component of the velocity field 5 . This allows us to develop the intensity and the velocity centroid maps, as 
well as their autocorrelation functions, according to the powers of cr p /po- It should be made clear that the 
expansion in Eq. i|13|) is not necessarily true for real data, but that it is used here as a reasonable first step 
towards understanding the effects of realistic density fluctuations on velocity centroids. 





Figure 2: Interpretation of the average M PtP appearing in Eq. (|23|l . The slab is viewed edge-on in the left 
figure, and face-on in the right. The value of M PtP at any lag R is the average of the correlation function B p p 
over all pairs of points whose projected separation on the plane of the sky is R. Three such pairs are presented 
on the left. 



5 It should be stressed that no hypothesis is made on the strength of velocity fluctuations, so that scaling them to the mean 
velocity vq is merely a convenient way to symmetrize expressions, and should not be given too much importance. We are aware, 
however, that the sound speed would be a much more physically meaningful velocity scale, and it will be used as such when the 
mean velocity is zero (see appendix [BJ. 
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4.2 Autocorrelation function of the intensity I 

Taking the autocorrelation function Aj(H) of the intensity map /(X), and using the linearity property of 
averages, we have terms of order up to two in density fluctuations, namely 

Aj(R) = (/(X)J(X + R)) = (a Po D) 2 [l + (y px ) + (y px+R ) + ( Mpx+R >] , (16) 

where we recall that the brackets stand for an average over X. These are performed on quantities integrated 
over z, and can therefore be interpreted as averages over the whole turbulent slab. Indeed, if one considers a 
three-dimensional field /(x), its integrated map F(X) defined by 

F(X) = 1 J /(X, z)dz has an average (F(X)) = ±J (/(X, z)) dz = ± JJJ /(x )dx (17) 

Id Id 

over a surface S of the sky. The quantity (F(X)) can then be seen as an average of / over the volume DS, and, 
assuming ergodicity, it can also be identifid with the ensemble average /(x) at any given position x within the 
flow. In the present case, the term (y Px ), which is of the first order in density fluctuations, can be written as 



PoD J po 

Id 



(18) 



since, by definition, the average of the density fluctuations over the whole volume of the slab is zero. The 
assumption of homogeneity allows one to write that (y pyi+R ) = as well. As for the second order term, which 
reads 



(2/pxJ/px+R.) 



J 6p(X, z)dz J Sp(X + R, z)dz \ = J J (Sp(K, «i)<5p(X + R, z 2 )) dz x dz 2 , 

Id J Lin J ' ° I 2 D 

(19) 

and where ify is the square domain Id x In, it includes an average over X that, following the idea used above, 
should be replaced by an ensemble average characteristic of the turbulent flow. Indeed, one can see that 



(Sp(X, zi)5p(X + R, z 2 )) = <5p( Xl )5p(x 2 ), (20) 

with the three-dimensional vectors xi = (X, z\) and x 2 = (X + R, z 2 ). Introducing the autocorrelation function 
B p p of the density fluctuations, which is defined, for any pair of points (x, x + r) within the slab, by 6 



B„,„(r)=Jp(x)5p(x + r), (21) 
we therefore can write the term under study as 

(yp x y PX+R ) = ^2 J J B pA x 2 - xi)dzidz 2 or, more concisely, (y px y px+R ) = ^M p , p (R), (22) 
il 

where M PiP (R), which is defined by the double integral 

M p , p (R) = -L J f j Bp . p (H + (z 2 - Zl )e z )d Zl dz 2 , (23) 



6 We chose to write B PiP instead of using the more consistent notation Ag p for the autocorrelation function of the density 
fluctuations, because, in the following developments, higher order and multipoint correlation functions of mixed fields arise for 
which the Ajr notation would have become cumbersome. 
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is interpreted as an average of -B PiP (r) taken over all pairs of points within the slab whose three-dimensional 
separation r has the given vector R for component in the plane of the sky (see Fig. |2J. Eventually, the 
autocorrelation function of the intensity simply reads 

A I (R) = (aD) 2 [p 2 +M p , p (R)]. (24) 

This expression will be exploited later on when compared with the autocorrelation functions of velocity centroids. 



4.3 Autocorrelation function of the non-normalized velocity centroid C 

The same general method applies when considering the non-normalized velocity centroid C(X). Its expansion, 
written in Eq. (|10(l . includes terms of order zero and one in density fluctuations. Therefore, its autocorrelation 
function Ac(R) features terms of order up to two. Namely, 

2 

A C (R) = (ap v D) 2 ([1 + y vx + y px + y pvx ] [l + y vx+R + y Px+R + y pvx+R ] ) = {ap v D) 2 ^ (a n ) , (25) 

71=0 

where a n is a term of order n in density. The explicit forms of these coefficients are given by 

fl l = VPX + ?/px+R + ypxVvx + R + y^xJ/px+R + VpVX + VpVX + R + ypVxUvx + R + VpVX + rVvx: (26) 
a 2 = UpxUpx + R I" VpVxUpx + R ' ypVx+RUpx ' UpVxVpVx + R- 

The linearity of the averaging process over the plane of the sky implies then that we should consider terms of 
the form (y\ x ) an d (2/AxZ/px+r)j where A and p represent either p, v or pv. Due to the homogeneity hypothesis, 
terms of the form (j/Ax+r) are of course equal to {yx x ), and so, considering the zeroth order contribution (ao), 
we have 

(a ) = 1 + 2 (y Vx ) + (yv x y-v x+R ) = 1 + — Sv(x) + —M VyV (R) = 1 + —M VtV (R), (27) 
with M VA ,(R) being defined similarly to M PjP (R) in Eq. (|23() . 

M ViV (R) = 1 1 B v , v (R+{z 2 - zi)e z )dzidz2, (28) 



and B VtV (r) = 5v(x)Sv(x + r) is the autocorrelation function of the fluctuations of the line-of-sight velocity 
component. Turning to the first order term, we have, since (y px ) = {y Px+R } = as shown in the previous 
section, 

< fl i) = (y Px yv X+R ) + (yv X y PX+R ) + 2 (y pvx ) + (y P v X y vx+R ) + (y P v X+R yv X ) (29) 

Each term in this equation can be related to correlation functions of the fluctuation fields 5p and 5v, beginning 
with 

(v Px yv x+R ) = ^t^j [[ ($P(X, 2i)5u(X + R, z 2 )) dz x dz 2 = ^5 // B p „(R + (z 2 - z 1 )e z )dz 1 dz 2 , (30) 

p v Q D^ JJ PoVoD 2 JJ 



introducing the mixed correlation function £? Pj „(r) = Sp(x)Sv(x + r). Considering the next term (?/ux?/px+r)j 
the homogeneity hypothesis allows us to shift the arguments by — R so that, exchanging the integration variables, 



(Vv x y PX+R ) = {yp x yv x - R ) = // B p>v (-R + {z x - z 2 )e z )dz 1 dz 2 , (31) 
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xI/dx+r) — n 2 [Bp,v(Ti- + (?2 zi)e z ) + B p<v ( R + (zi Z2)e z )] dzidz2- (32) 



Combination of both terms leads to 

p u D 2 

Now, according to iMonin fc Yagloml l)l975h . and assuming that the four-dimensional field (p, Sv), made up of 
the density p and vector velocity Sv, is homogeneous and isotropic, the correlation function of p and 6v is a 
vector quantity, and due to isotropy should be along the separation vector r. It should therefore be of the form 
p(x)5v(x + r) = f(r)r, so that, projecting this relation on e z , we have 



f{r)r.e z = p(x)5v(x + r).e z = <5p(x)<5v(x + r).e z = fy(x)5u(x + r) = S p ,„(r), (33) 



using the fact that <5v(x) =0. As a consequence, we have B P}V (— r) — —f(r)r.e z — —B PjV (r), so that B p _ v is 
an antisymmetric field. Of course, this is a statistical property, and does not hold when considering the specific 
values of velocity and density fluctuations for a given pair of points within the flow. What one can conclude is 
that the integrand in Eq. (|32|) is zero, and so {y vx yp x+R ) + (2/px^ux+r) = 0- So is the next term (y pvyL ), since 

(Vpv x ) = I (W, z)Sv(X, z)) dz = - J— f B p . v {Q)dz = 0, (34) 

PqVqD J Pqv D J 

Id Id 

because the antisymmetry of B p v implies that B PiV (0) — 0. Similarly, the last two terms of (a±) are given by 

x+rZ/^x) 2n2 // [-^pVtVyR + {z2 - Zi)e z ) + B pv>v (-R + (zi - z 2 )e z )) dzidz 2 , (35) 

PoV -L> JJ 
T 2 

1 D 



introducing the two-point correlation function B pv _ v (r) = 5p(x)Sv(x)Sv(x + r) . This last combinat ion of terms, 
unlike the one considered previously, is not necessarily zero. According to lMonin fc Yagloml l)l975|) . it is of the 
form 

B pVtV (v)=g(r)^^- + h(r), (36) 



where g and h arc functions of the scalar separation r. This form stems from the fact that Sp(x)Svi (x)Svj (x + r) , 
where i and j stand for any of the x, y and z coordinates, is a tensor of rank 2 which we suppose to be 
homogeneous and isotropic. It follows that B pv ^ v is symmetric, and so, using an averaging notation M pVtV 
defined, similarly to the expressions of M PtP and M v v , by 



M pv , v (R) = ^ II B pv ^(R+(z 2 - z 1 )e z )dz l dz 2 , (37) 



n, 

we can write (ai), contribution of the first order in density fluctuations to the autocorrelation function Ac, as 

1 2 

(ai) = (y P v x yv x+R ) + (y P v x+R yv x ) = — 2 [M P v,v(R) + M pVyV (-R)} = — ^M pv , v (R). (38) 

Po v o Pov 

Lastly, the second order term in density fluctuations (02) = (y Px y Px+Ft .) + (v 'pv x V Vx+r) + \V 'p^x+nV 'px) + 
(ypvx.ypv x+R ) is computed in much the same way, introducing the two-point correlation functions 



B pv , p (r) = Sp(x)Sv(x)Sp(x + r) and B pViPV (r) = 8p(x)8v(x)6p(x + r)Sv(x + r), (39) 
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and their averages defined by 
1 



M, 



p«,p(R) — jj 2 



B pVtP {R + (z 2 - zi)e z )dz 1 dz 2 and M pVtPV (R) = — 



B pv<pv (R + (z 2 - z 1 )e z )dz 1 dz 2 . 

(40) 



With these notations, it is straightforward to derive the expression of the second order term, 



(02) = -2 



M P}P (R) 



^M$ p (R) + ±M pv , pv (K) 



(41) 
(42) 

Now, still following iMonin fc Yagloml l)l975(l . B pv ^ p exhibits the same antisymmetry property as B PiV and we 
can conclude that the symmetric part of M pViP (R) is zero, so that the autocorrelation function Ac of the 
non-normalized velocity centroid eventually reads 



with the (s) superscript indicating the symmetric part of a given function, defined by the relation 

/«(!) = J [f(x) + /(-*)]. 



A C (R) = (aD) 2 [p 2 v 2 + p 2 M v , v (R) + 2p M pv , v {R) + v%M p , p (R) + M pVyPV {R) 



(43) 



When considering only the zeroth order in density fluctuations, one finds that the expression of Ac(R) is very 
similar to that of Aj(R), with Ac(R) — (aDpo) 2 [v 2 +Af U) „(R)l. This is not surprising, as we shall discuss 
in section [S] Moreover, as the forthcoming comparison of the autocorrelation functions of both normalized and 
non-normalized velocity centroids will be performed on expressions of order up to one only, it is useful to write 
out the truncation of Ac(R) at that order, which is, given Eq. I|43|l . 



A C (R) = {aDf [p 2 vl + p 2 M VyV (R) + 2p M p ^(R)] . 
4.4 Autocorrelation function of the normalized velocity centroid Co 



(44) 



There is little qualitative change to the method when dealing with the normalized centroid. Contributions of 
each order can be computed as easily as for the non-normalized case, with the main difference being that, given 
the expansion written in Eq. (|13|l . Ac {R) is a theoretically infinite series which reads 



-o 2 



m n— 



with bri 



p=0 



, (45) 

the generic term (a n b m ) obviously being of order n + m in density fluctuations. Now, given that the highest 
order present in Ac(R) is two, it seems reasonable to consider only terms of order at most two in the expansion. 
However, it may also be argued that, since density fluctuations effectively contribute to the first order in 
the expression of Ac{R), assessing the effects of normalization may already be performed when limiting the 
expansion to that order, 

A Co (R) = vl [(a b ) + (aoh) + (aib )] = v 2 (S + Si) with S Q = (a b ) and Si = (aoh) + (aib ) (46) 
Now, the zeroth order term in Ac (R) has been computed in the previous subsection, 



So = (a b ) = (a ) = 1 + —M ViV (R), 



(47) 
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so we turn to the first order term S± = (ao&i) + (ai&o) which reads, after a long but straightforward calculation, 



Si = 2 (y pvx ) - 2 (y Px y Vx ) + (y P v x yv x+R ) + (y )-(y >• (48) 

The first term is zero, as shown before. So is the second term, for we have 

(y P xVv x ) = ^2 f[ B p,v{{z2 - zi)e 2 )dzidz 2 = 0, (49) 

PoVqU J J 

1 D 

since the integration domain is symmetric about z 2 — z\ — and B p ^ v is antisymmetric. The combination of 
third and fourth terms has been computed in the previous subsection 

2 

(j/pfx^x+R,) + (2//tox+r2/«x) = 2 -Mpi;,tj(R)- (*^) 

P0 V 

The remaining terms require a more elaborate, although similar, treatment. Taking for instance the last term, 
we have 

(yv x yv x+R y Px+R ) = 2^3 JJJ (<M X ) zi)5v(X + R, z 2 )Sp(X + R, z 3 )) dzidz 2 dz 3 , (51) 
° i 3 

1 D 

where the integration is performed over the cubic domain I D = Id X Id X Id- Introducing the three-point 
correlation function B VyVp (ri, r 2 ) = dv(x)5v(x + Ti)Sp(x + r 2 ), this expression reads 

(yv x yv x+R y Px+R ) = tt^ [[[ B v . v . p (R + (z 2 - Zi)e z ,R + (z 3 — z 1 )e z )dz 1 dz 2 dz 3 = —^M v . VyP (R, R), (52) 

PoVq^ JJJ Po v o 

1 D 

with an averaging notation M v<v P (R, R) reminiscent of the ones used previously, 

M v<ViP (R, R) = i JJJ B ViV>p (IL + (z 2 - Zi)e Z! R + {z 3 - z 1 )e z )dz 1 dz 2 dz 3 . (53) 
i 3 

1 D 

Similarly, the next-to-last term can be written as 

1 2 

(2/TOC^x+R.ypx) = 2 Af *Wp(- R ) - R ) S0 that (^x^x+rJ/px+r) + (^x^x+r^Px) = 2 M S.p( R ' R )> 

PO^o P0 U 

(54) 

with Mif2,p(R, R) being the symmetric part of the function M v<v p (R, R), 

M« >p (R, R) = i [M Uit ,, p (R, R) + M„,„ ;P (-R, -R)] . (55) 

Finally, the autocorrelation function Ac (R) of the normalized centroid has the following expression when 
limited to contributions of the first order in density fluctuations, 

A Co (R) = vl + M v , v (R) + — \M pv , v (R) - M« (R, R)l . (56) 

Po L J 

It appears then that normalization of velocity centroids indeed performs a first order correction, as compared 
with Eq. 144fl . although it does not necessarily fully remove the density structure, as will be discussed in more 
detail in section 
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5 The case of fractional Brownian motion fields 



Explicitly computing the forms of the two-dimensional statistical measures Ai (R) , Ac (R) and Ac (R) requires 
the knowledge of the three-dimensional correlation functions appearing as integrated terms in equations (|24[1 . 
()44|) and l|56p. However, it may prove impossibly difficult to build a complete and consistent set of such functions 
and then to analytically compute their averages. Therefore, as a first step, one should stick to the lowest order 
terms in the expressions of the two-dimensional autocorrelation functions, 

A I (R) = (aD) 2 [p 2 + M p , p (R)], A C (R) ~ (aD Po ) 2 [v 2 + M VjV (R)] and A Co (R) ~ v 2 + M„,„(R). 

(57) 

In this limit, one only has to compute the M p>p (R) and M VyV (IL) averages given density and velocity fields 
with known autocorrelation functions B PtP (r) and B VjV {y). One possibility is to suppose that both of the three- 
dimensional fields are fractional Brownian motions (fBm) . These are defined in the following way: if T is such 
a field in N dimensions, it is characterized by a single-index power-law structure function, 



5^(r) = [^(x + r)-^(x)] 2 = 2A(-j , (58) 

where r = |r| is the length of the separation vector, A is a positive constant and H is a real number in [0, 1] called 
the Hurst exponent. The restriction H ^ corresponds to the fact that the amplitude of fluctuations should 
decrease at smaller scale, while we impose H 1 because an fBm field having a Hurst exponent H > 1 would 
be uniform, as shown i n app endix IdI F ractional Brownian m otion fields have already been used l|Stutzki et all 
1998| iBrunt fc Heverl . l20f)2t iMiville-Deschenes et all l2003al) to model clouds of the diffuse, non-starforming 
interstellar medium. In this section, we shall use them to model the three-dimensional density and velocity 
fields, so that, for the former, 

S p (r) = [p(x + r)-p(x)] 2 = [<V(x + r)-,Sp(x)] 2 = 2 [a 2 p - B PtP (r)] = 2A P (£) , (59) 

since the density fluctuations field is supposed to be homogeneous. Similarly, for the velocity field, 

S v (r) = [«(x + r)-«(x)] 2 = [Sv( X + r)-Sv( X )} 2 = 2 [a 2 v - B v ^{y)] = 2A V (^)^ ■ (60) 

These last two equations make use of the relationship between the second order structure function of an 
homogeneous field T and its autocorrelation function Ayr, which is Sjf(R) = 2 [A^-(O) — A^R)]. In the results 
of sectional we then have to inject the following relations 

and B VlV (T)= *-A v (p) ■ (61) 

To compute the averages M„ „( R) and M n n (H) featured in the e xpressions of ^4/(R), Ac(R) and Ac (R), we 
can use the calculation scheme of lChaudrasekhar fc Miiur hl ill 952^ . which is given in appendix 1X1 to write them 
as single integrals over the separation Az — z-i — z\ along the line of sight, noting that Az G I2D = [—D,D]. 
We then have 

M VtV (R) = -^2 J(D-\ Az\)B VtV (R+Aze z )dAz and M PtP (R) = J(D-\Az\)B PtP (K+Aze z )dAz. (62) 

In these expressions, the factor (D — \Az\) represents the weight of each separation Az, that is the relative 
number of pairs of points taken into account whose projected separation along the line of sight is Az. With the 



r 



2H 
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expressions above for B v ^ v and B P:P: it is possible, as shown in appendix El to write these averages under the 
form 

M ViV (H) = -A v K(R,H v )+a 2 v and M PjP (K) = -A p K{R, H p ) + a 2 , (63) 
where R = |R| is the scalar separation on the plane of the sky, and the function K is given by 



K(R,H) = 



D 



(R 2 



z dz 



1 



H + l 



R*_ 

D 2 



2\ H+l 



D 2 



2\ H+X 



(64) 



The integral in Eq. I|64|) cannot be explicited ijGradshtevn fc RvzhikL Il980|) unless H — or H = 1, but one 
interesting limit to consider is that of small separations on the plane of the sky (R -^i D), in which case K(R, H) 
can be developed in powers of R/D. As shown in appendix IU1 this yields the following scaling relations for 
M VtV (R) 



M„,„(R) aai- A v a(H v ) - A v b(H v ) — 



M„,„(R) ~ ai - A v a{H v ) - A v b{H v ) — 



M v , v (B.)^a 2 -^+A v (^) ln(- 



R 


v 2H 


D 








( 


i) 


In 











for < H v < -, 



for - < H v < 1, 



for H v = -, 



(65) 
(66) 
(67) 



where a and b are functions of H v only. Similar relations are satisfied by M p p (R) depending on the value of 
H p . To the lowest order, the structure function 5c (R) of the normalized centroid Co therefore reads 



S Co (R) = 2A v b(H v ) 



S Co (R) = 2A v b(H v ) 



2H v + l 



for 0^H V < 



2" 



for - < H v < 1, 



^c (R) 



R 



D 



2A„ — In - 



1? 



for H v = -. 



(68) 

(69) 
(70) 



According to Eq. I|57|l . the structure function Sc of the non- normalized velocity centroid has the same scaling 
behaviour. This is also the case for the structure function Si of the intensity, depending on the value of H p . 
The consequences of such forms are presented in the next section. 



6 Discussion 

As expected, the statistical measures on the intensity and centroid maps can be expressed in terms of the 
statistical properties of the underlying three-dimensional velocity and density fields. More precisely, it should 
come as no surprise that the autocorrelation functions Aj (R) , Ac (R) and Aq (R) should be written as averages 
of correlation functions within the slab in the manner described in Fig. [21 since all pairs of points with a given 
separation R in the plane of the sky should contribute to the two-dimensional statistical measures at lag R. 
It also appears that the zeroth order term in the autocorrelation functions of the velocity centroids has the 
same form as the autocorrelation of the intensity map, since we then have 

A C (R) ~ (a Po D) 2 A Co (R) ~ (ap D) 2 [v 2 + M^R)] while A/(R) = (aD) 2 [p 2 + M p>p (R)] . (71) 
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The identity of these forms was to be expected, since the zeroth order terms in the autocorrelation functions of 
the velocity centroids are the limits obtained when the density is uniform. In this case, the centroids simply are 
integrals of v, as the intensity is a simple integral of p. In the context of fBm fields, the implication of this limit 
is that both non-normalized and normalized velocity centroids have a fractional Brownian motion behaviour, 
since their structure functions are power laws, with respective Hurst exponents He and Hq such that 



He — Hen — H„ 



- for < H v < - 
2 2 



and H c = H Co = 1 for \^<H V 



?S 1. 



(72) 



And similarly, the intensity map has a fractional Brownian motion behaviour with a Hurst exponent Hi with 



Hi — H n 



— for 
2 



< H p < \ 



and Hr = 1 for - < H a < 1. 

2 p 



(73) 



These results should be interpreted in the light of what is already known about the statistical properties of the 
integrated emission map / from studies in the Fourier domain. Indeed, in the case of optically thin lines, the 
po wer spectrum ind ex 7/ of the intensity map is the same as that of the three-dimensional density field 7 P (see 
e.g. iGoldmanl l2000|) provided the depth probed D is larger than the transverse scales R, which is the case in 
the limit considered in section [3J Now, spectral index 7 and Hurst exponent H are related by 7 = 2H + N, 
with N the dimension of the space over which the field is defined. As a result, IStutzki et alJ l)l998j) concluded 
that the Hurst exponent of the integrated emission map is Hi = H p + 1/2 for H p ^ 1/2 with a turnover to 
ffj = 1 for H p ^ 1/2, since 



Hr 



II 



lp 



H n 



1 

2' 



(74) 



using the fact that the integrated emission map is two-dimensional, while the original density field is three- 
dimensional. This result is precisely what we find from our analytical study performed solely in real space. 
Simi larly, the expression for the st ructure function of the normalized centroid confirms the numerical findings 
of iMiville-Deschenes et afl l)2003aj) . who showed that the spectral index 7c of normalized centroid maps was 
equal to that of the three-dimensional velocity field, j v . It should be noted that their result holds even for 
density fields with large fluctuations, which implies that our analytical result may be applicable to a more 
general class of fields. On the other hand, in the limit of neg ligible d ensity fluctuation s, one should find the 
same behaviour for the non- normalized centroid, a result :Miville-D eschenes et al have not reported on. 

Returning to a study of the autocorrelation functions of velocity centroids in the first order approximation, 



AeCR) 



■ Af„,„(R) + — M p ^(R) 



Po 



and A Co (R) 



■M V<V (R)- 



2 



(R) 



M«,,(R,R) 



we see that, since M pv ^ v (R.) — M^2 !P (R, R) is not necessarily zero, normalization actually does not remove 
the first order contribution of density fluctuations. The empirical assumption that normalization somehow 
eliminates the influence of density fluctuations on velocity centroids may therefore not be true, although assessing 
the magnitude of the remaining first order contribution in the normalized centroid's autocorrelation function 
with respect to the non- normalized case could be difficult, as it involves computing averages of correlation 
functions whose forms are still unknown. This will be investigated in a forthcoming paper. 
Such shortcomings of the normalized centroids as a way to retrieve the act ual velocity statistics from the 
observational data have already been pointed out bv lLazarian k, Esauivell l)2003|) . albeit in the case of simulations 
of highly compressible MHD turbulence. Instead, they introduced modified velocity centroids (MVC) C m , 
defined through their second order structure function Sc m as 



5 Cm (R) = ( [C m (X + R) - C m (X)] z ) = 5 C (R) - K + < + <)Si(R 



(75) 
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where Sc and Sj are the structure functions of the non-normalized velocity centroid and of the intensity, 
respectively. Given the relationship between the structure and autocorrelation functions, we can write the 
autocorrelation function Ac m of the modified velocity centroids as 

A Cm (R) = A C (R) - (v 2 + a 2 v + a t 2 h )^/(R) + E, (76) 

E being a constant involving the values of Ac and Aj at zero lag. Now, since Aj contains no first-order 
term, it appears that modified and non-normalized velocity centroids have the same first order contribution 
in their autocorrelation functions. As noted earlier, however, it is not yet clear whether this contribution is 
actually more important than the first order term in the autocorrelation function of normalized centroids. If 
it is, it would tend to prove that normalization may be a better way of retrieving velocity statistics in weakly 
compressible flows. To make the comparison of correction schemes clearer, one can show a relationship between 
the autocorrelation functions of /, C and Co as written in equations 124|) . I|43|l and (|56[1 . 

{aD Po ) 2 A Co (R) = (aDp vo) 2 + A C (TL) - «gi4/(R) - 2(aD) 2 ^M« iP (R, R) + G(R), (77) 

where G(R) stands for the terms of order at least two in density fluctuations. Deriving the relationship between 
the structure functions of /, C and Co, the latter being noted Sc , we have 

(aDp Q ) 2 Sc (R) = Sc(R) - ugS/(R) - ^D) 2 Po [m«, p (0, 0) - M«, P (R, R)] + 2 [G(0) - G(R)] . (78) 

Comparing this with Eq. I|75|l , it is clear that modified velocity centroids subtract more of the intensity structure, 
thus removing second order terms, while normalization performs a more complex correction, involving all orders 
in density fluctuations. Actual comparison of the me rits of normalized and modi fied velocity centroids requires 
numerical tests, which have only been performed bv lLazarian & Esauivell l^flO.'j ) on highly compressible MHD 
turbulence. That particular case is actually not within the scope of our study, firstly because density fluctuations 
are large and the condition for the expansion in the normalized velocity centroid's case may not be met, and 
secondly becaus e the presence of a magnetic field implies anisotropy. Both limitations render the comparison 
with the work of lLazarian fc Esauivell l|20Q3) a bit uncertain, although it is likely that terms of higher order in 
density fluctuations may become dominant. Thus, the effectiveness of modified velocity centroids, as compared 
to normalized centroids, in their simulations may be related to the second order correction through Si, which is 
more important in modified velocity centroids. In any case, further analytical and numerical work is warranted 
in order to establish these conclusions more firmly, especially in the weakly compressible flow regime. 



7 Conclusions 

An analytical study of a simple slightly compressible turbulent cloud model was presented, assuming homogene- 
ity and isotropy of the turbulent flow. From the expressions of the antenna temperature for an optically thin 
spectral line and of its successive moments with respect to the line of sight velocity component, we computed 
the autocorrelation functions of the intensity and of both normalized and non-normalized velocity centroids, 
which involve averages, along the line of sight, of correlation functions of the three-dimensional density and 
velocity fields. 

To the lowest order, the autocorrelation functions of the velocity centroids behave, with respect to the velocity 
field, as the autoco rrelation function of the i ntensit y with respect to the density field. This sheds light on the 
numerical result of iMiville-Deschenes et all 1 2003al) . who found that, for fractional Brownian motion density 



and velocity fields, the spectral index of the normalized centroid is equal to that of the velocity field. We derived 
this result analytically, for separations across the sky much smaller than the cloud's depth, and in real space, 
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w hile previous studies such as tha t of iGoldmarJ l)2000|) were performed in Fourier space. However, the result 
of lMiville -Deschenes et al.l l)2003al) holds for fields outside of the validity domain for our calculation. 
Comparison of the expansions of the autocorrelation functions of both types of velocity centroids shows that 
normalization performs a correction of the first order in density fluctuations, although its magnitude remains 
to be assessed, a task for which numerical simulations are probably necessary. Numerical tests should also 
pr ovide us with a robus t comp arison between normalized velocity centroids and the modified velocity centroids 
of lLazarian fc Esauivell |2003), which imply corrections of order two in density fluctuations. At present, this 
comparison has been performed only on simulations of highly compressible and magnetized turbulence, a case 
beyond the scope of our analytical study, and has shown that, in this particular case, modified velocity centroids 
provide a more reliable tool than normalized centroids. 

In a forthcoming paper, we shall therefore present numerical simulations aimed at assessing the validity domain 
of our calculations and, beyond normalized and modified velocity centroids, pursuing the search for a better 
correction scheme able to retrieve the underlying velocity statistics from observational data. 
I wish to acknowledge fruitful discussions with Alex Lazarian during his stay at the Ecole normale superieure. 
His suggestions, in the early stages of this work, proved very helpful. I also wish to thank Enrique Vazquez- 
Semadeni for his careful reading of the manuscript and insightful remarks. 
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A The Chandrasekhar-Miinch scheme 

For a function /(r) depending on the three-dimensional separation r, double integrals of the form 



J= II f(R+(z 2 -z 1 )e z )dz 1 dz2 (79) 

can be transformed, taking into account the fact that the integra nd only depends on the difference Az — Z2 — Zi, 
because R is fixed, as shown in C handrasekhar fc Mu nch ( 1952). Making the variable change Az = Z2 — z\ and 

Z = Zi, 

1 = JdzjdAzf(R + Aze z ). (80) 



In In, 



where the notation lr> z stands for the segment Id shifted by —z, Io z = [-D/2 — z,D/2 — z], which we can 
separate in two, with a negative part, 1^ = [-D/2 — z, 0], and a positive part 1^ = [0, D/2 — z\. 

1= fdzj dAz /(R + Aze z ) + JdzJ dAz/(R+ Aze z ). (81) 
Id ^ Id IB, 

Exchanging the order of integrations in each of these two terms, and expliciting the integral over z, we have 

1= f(H + Aze z )(D + Az)dAz+ / /(R + Aze z )(D — Az)dAz, (82) 

J-D JO 
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which can be finally written as 



1= J(D -|Az|)/(R + Aze z )dAz. 



(83) 



The factor (D — |Az|) expresses the fact that for a finite slab, two given separations Az and Az' are not 
represented by the same number of pairs of points. 

B Autocorrelation functions of the velocity centroid maps in the 

case vq = 

When the mean velocity vq is null, the method used in the main body of the paper cannot be applied without 
some modification. Little needs to be done, however, as velocity fluctuations can be measured with respect to 
the sound speed c s , which we assume to be uniform within the turbulent slab. We may then write the velocity 
centroid C as 



C(X) = aj(p a + 5p(X, z)) Sv(X 7 z)dz — ap c s D 

Id 



1_ f Sv(X,z) 1 f Sp(X,z) Sv(X, z) 



Id 



Id 



Pa 



(84) 



which can be written in the shorthand form C(X) = apoc s D [y Vx + Vpv x ]- The normalized velocity centroid 
then reads 

Co(X) = ^ = c s [y vx + y pvx ] [1 + y^ 1 . (85) 

The calculation therefore proceeds in the same way, with c s taking the place of vq and the of Eq. (|26|l being 
replaced by new coefficients a[ given by 



a o — yox!/«x+Ri 

a l — 2/p-ux2/i>x+R "I" J/phx+rJ/dx ) 
a 2 = ypfx2/pi>x+R- 



It follows that the autocorrelation function of the non-normalized centroid is, in this case, 

A C (R) = (aD) 2 [p%M VtV (R) + 2p M pv , v (R) + M pv , pv (R)] , 



(86) 



which is precisely what is found form the general case when setting vq — 0. Similarly, for the normalized 
centroid, 



A Co (R) = M V>V (R) + — M P „,„(R) - M^ s 2 „(R, R) 
Pa L 



(87) 



C Expansion of the averaged autocorrelation function of an fBm 
field 



We wish to compute the following expression as a function of the lag R = |R|, in the limit D, 



M(R) 



1 

W 2 



-A- 



R A + z z 



D 



2H 



(D-\z\)dz, 



(88) 
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where A is a positive constant and H € [0, 1] is the Hurst exponent. M(R) can be expressed as an integral over 
[0,D], 



M(R) 



D 2 J 



-A 



(i? 2 + z 2 ) 



H 



D 



2H 



(D - z)dz 



D 2 



-A 



D (R 2 + z 2 



H 



-(D-z)dz + a z / (D-z)dz 



The last integral is easily shown to be equal to D 2 /2, and so 



M(R) = - 



2A 



JJ2H+1 



(R 2 + z 2 ) H dz 



2A 



£)2H+2 



[ L \r 2 + z 2 f zdz + o 2 
Jo 



(89) 



The second of the remaining integrals can be explicited by setting R 2 + z 2 = r 2 , since zdz = rdr for R fixed, 



(R 2 + z 2 ) H zdz 



1 



s/R 2 +D 2 

r 2H+1 dr = 
r + 2 



(R 2 + D 



2 \H+1 



(R 2 ) 



2\H+1 



since 2H + 1 ^ — 1 for ^ H ^ 1. We therefore have 

rD A 



H + l 



l + ~D 2 



D 2 



2\ 



For a zero separation in the plane of the sky, this expression is, explicitly, 



M(0) = 



2A 



JJ2H+1 



z 2H dz 



A 



H + l 



a 2 = - 



A 



(2£r + i)(ir + i) 



(90) 



(91) 



(92) 



When R ^ 0, the first integral in Eq. (|91|l generally cannot be written in closed form ijGradshtevn fc RvzhikL 
Il980l) . Notable exceptions are for H = and H = 1, when, respectively, 



M(R) = -A + a 2 and M(R) = -A 



A 



(93) 



To simplify the notations, we introduce the function K = K' — K", defined by M (R) = —AK(R, H) + a 2 , with 



K'(R,H) = ^L T J\r 2 + z 2 ) H dz and K"{R,B) = -±^ 



If 

D' 2 



2\ H+l 



2\ H+l' 



(94) 



It is now interesting to consider the case of separations R -C D, which corresponds to studying the small 
scale structure of centroid maps. This allows to develop K(R, H) in powers of R/D. Expansion of K" is 
straightforward, 



K"(R, H) = —^—Y 



J n (H + l) (R 



n>0 



D 



2/i 



for 0<H^1 and K"(R,0) = 1, 



(95) 



where we introduced J n ( x ) = x{x — 1) . . .(x — n+1). As for A"', since we consider D > i? > we can write it 
as 



fD/R 

K'(R,H) = 2 (l + y 2 ) H 

Jo 



= 2 



2ff+l 



D/R 



H 



(96) 
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The first integral, between and 1, which we dub Kq(H), cannot be explicited, but since it does not depend on 
R, it is unimportant with respect to the structure of the moment maps. The second integral can be transformed 
in order to develop the integrand through l + y 2 = y 2 (l + y~ 2 ), 



h"[R.II) >(| 



2H+1 



(97) 



Computing the integrals in the equation above for < H < 1 and H ^ 0.5, we have 



K'di.ij)- 2(| 



2H+1 



k>0 



2H-2k+l 



(98) 



which, after rearranging the terms, leads to the following expression for K' , 

2H+1 



K'{R,H) = 2 



Ko(H) 



lk(H) 



fe^O 



(2H - 2k + l)k\ 



E 

k>0 



R 



(2H-2k + l)k\ \D 



2k 



(99) 



If < H < 0.5, then 1 < 2H + 1 < 2 and the leading order expansion is therefore 



K'(R,H) 



2H+1 



+ 2 



lk{H) 



2H - 2k + 1 



2H+1 



(100) 



while if 0.5 < H < 1, then 2 < 2H + 1 < 3 and the k — 1 term of the last sum is dominant, 



K'(R,H) 



2H ( R 



2H + 1 2H- 1 \D 

This last expression is also valid for H — 1, while obviously K'(R, 0) = 2. For H = 0.5, the k = 1 integral is 



(101) 



*o(0.5)-± 



and for k ^ 1 the integrals are unchanged. It follows that 

K'(R,0.5) = 2 
which gives, after rearranging the terms, 

R -' (R ,„.5> = l-(^) 2 ,n(£) + 



7fc(0-5) 
2-2fc 



2-2k 



(102) 



(103) 



i? 2 

D 2 



\ - 7fc(0-5) (R 
fri 1-k \D 



k>2 



2k 



(104) 



Finally, the leading order expansion of K(R, H) for small separations i? <§C D is 

2JJ+1 



K(R,H) 



(2JT +!)(# + !) 



+ 2 



^-Es^^m (I) '» °^<°- 5 - 

fe>0 v 7 



(105) 
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K(R,H) 



1 ( R 



(2H +1)(H + 1) 2H-l\D 



for 0.5<i/^ 1, 



(106) 
(107) 



Consequently, the expansion of Af(R) = —KK{R 1 H) + a 2 reads 



M(R) ~ cr 2 



(2JT+1)(# + 1) 
M(R) ~ a 2 - 



2A 



fe^O 



2H-2k+l 



I) 



2H+1 



for O^H < 0.5, (108) 



A 



(2if+ 1)(# + 1) 2£T-1 V£> 



A f R\ 2 

I — I for 0.5<if^l, 



M(R)~<^--+A 



(109) 



(110) 



The consequences on the statistical measures performed on the moment maps are straightforward to derive. 
They are given and analyzed in the main body of the paper. 



D The uniformity of fBm fields with H > 1 

We consider the mean squared increment of an fBm field T of Hurst exponent H > 1 between positions x and 
x + r. The separation vector r is then separated in p equal parts, so that 



[^(x + r)-^(x)] 2 = ]T 
U=o L 



T I x H r 

P 



T ( x + -r 

P . 



p-i 
Lfe=0 



(111) 



where AF k is the increment of T between the k th and (k + l) th positions. Expansion of the expression above 
yields 

p-i _ 

E ( 112 ) 

k=0 i<j 



[^(x + r) - ^(x)] 2 = ^1 + 2 E A^A^, 



Now, ATiATj can be written as an autocorrelation product of the function Gr, P = J 7 (x + r/p) — J 7 (x), 



Gr,p ( x + -r 



Jr.p 



x H — r 



Since autocorrelation functions are decreasing from their zero spacing value, we have 



A^A^- ^ [G riP (x)] z , 
which allows us to give an upper limit for the mean squared increment of T , 



2A 



2A 



[^(x + r) - ^(x)] 2 < ^ l r 2H + p(p - l)^ 211 = 2kp 2i - 1 - H \ 2H . 



(113) 



(114) 



(115) 



This result being valid for any value of p, which characterizes a subdivision of the separation vector r, we see 
that since H > 1, the limit p — > oo implies [.F(x + r) — ^(x)] 2 = and therefore T is a uniform field. 
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